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Abstract 

The peak effect in weakly pinned superconductors is accompanied by 
metastable vortex states. Each metastable vortex configuration is charac- 
terized by a different critical current density Jc, which mainly depends on 
the past thermomagnetic history of the superconductor. A recent model [G. 
Ravikumar et al, Phys. Rev. B 61, R6479 (2000)] proposed to explain the 
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history dependent Jc postulates a stable state of vortex lattice with a critical 
current density J**, determined uniquely by the field and temperature. In 
this paper, we present evidence for the existence of the stable state of the 
vortex lattice in the peak effect region of 2H — NbSe2- It is shown that this 
stable state can be reached from any metastable vortex state by cycling the 
applied field by a small amplitude. The minor magnetization loops obtained 
by repeated field cycling allow us to determine the pinning and "equilibrium" 
properties of the stable state of the vortex lattice at a given field and tem- 
perature unambiguously. The data imply the occurence of a first order phase 
transition from an ordered phase to a disordered vortex phase across the peak 
effect. 

PACS numbers :64.70.Dv, 74.25.Ha, 74.60.Ge,74.60.Jg 
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I. INTRODUCTION 



In the presence of strong pinning, the vortex state of type II superconductors is usually 
characterized by the critical current density Jc{H, T) which decreases monotonically with 
increasing field H or temperature T. In the weakly pinned superconductors, on the other 
hand, the interplay between the intervortex interaction and the flux pinning produces an 
anomalous peak in J^, as a function of both field and temperature just below the normal 
state boundary (usually designated as the peak effect or PE). Within the collective pinning 
description 0, this signifies that the vortex phase undergoes a transition/crossover from 
an ordered state to a disordered state The detailed nature of this transition, e.g., 

whether it is a thermodynamic phase transition or not, remains a subject of considerable 
debate. 

One of the key issues is the detection of an anomaly in the thermodynamic quantities, 
such as, specific heat or equilibrium magnetization Meq. Jc and Meg can be estimated from 
the measured irreversible magnetization data of a superconducting sample f^-^ using the 
relations. 



MH) = [M{H i) - M{H ^)]/2gftoR, (lo) 

Meg{H) = [M{H t) + M{H i)]/2, (16) 

where M{H | ) and M{H |) are the magnetization in the increasing (forward) and decreasing 
(reverse) field cycles respectively, /iq = 47r x 10^^ W/A.m, R is the sample dimension 
transverse to the applied field and g is a factor which depends on the sample geometry. 
Eq. 1 implicitly assumes that Jc is history independent and is thus uniquely determined 
by the local induction B. However, across the peak effect region, the above equations 
are not valid due to a strong history dependence in Jc P,p|-p^. Recently, considerable 



efforts have gone into ascertaining the equilibrium magnetization across the peak effect 
region, where an order- disorder transition occurs in the vortex matter. However, these efforts 
have met with ambiguous and somewhat conflicting results. For example, the construction 



of the equilibrium magnetization from the hysteresis loop by using two different kinds of 



minor magnetization curves |]I6|,|T8],|T9| , results in apparently different conclusions. In one 



case, a jump ||16|,|19| in Meg could be found at the onset of the PE, while the other case 



shows no increase at all ||T8[. These differences apparently originate from the difficulties 
in establishing an unambiguous and reproducible vortex state due to a strongly history 
dependent configuration of the vortex matter in the PE region. The different procedures 
proposed to obtain Meg shall be discussed in section II. 



In Sec. Ill, we briefly discuss a recent phenomenological model |2^, which addresses 
the issue of the history dependent Jc and the metastability in the vortex state through an 
extension of the Bean's critical state model H. In Sec. IV, we present an experimental 



method based on the ideas of the model to obtain a unique "stable" vortex state in the 



PE region, which is independent of the past magnetic history. We propose that this state, in 
effect, is the "stable" or "equilibrium" state and evaluate the critical current density and Mgg 
of this state. We further demonstrate that, a sharp change in the equilibrium magnetization 
(albeit smeared) occurs across the PE region. These results imply that an underlying first 
order phase transition, presumably driven by a competition between elastic and pinning 
energies in a situation where thermal fluctuations are weak, marks the peak effect. 

II. MINOR CURVES AND THE EQUILIBRIUM MAGNETIZATION ACROSS 

THE PEAK EFFECT 

In the peak effect region, the critical current density in the increasing field cycle Jc{H f) 
is less than that {Jc{H [)) in the decreasing field cycle [p|, p!0| , p!4| for H < Hp, where Hp is the 



field where Jc is maximum. However, well below the onset of the PE and aX H > Hp, Jc is 
independent of the magnetic history. One of its consequences is the peculiar behavior of the 
minor magnetization curves, which can not be reconciled within the critical state model . 
For instance, a typical minor magnetization curve (type I) initiated from a field H < Hp in 
the PE region saturates without meeting the reverse magnetization curve [|H|,|IB|,|T7[ , as shown 
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in Fig. 1(a). On the other hand, the minor curves (type II) measured by increasing the 
field from different points on the reverse magnetization curve overshoot the forward curve 



1^,|T8[ as shown in Fig. 1(b). The two types of anomalous behavior may be contrasted with 



the conventional behavior for the minor curves starting at (a) H > Hp and (b) H<<Hpi, 
i.e., for fields well below the PE region. The latter catagories of minor curves meet the 
magnetization envelope, constituted by the forward and reverse curves, as expected from 
the Bean's critical state model. 



A new procedure was proposed by Roy and Chaddah |16| to obtain M^g from the minor 
magnetization curves of the type I by the relation. 



M,g{H) = [M{H + 6,]) + Mml{H - 6, i)]/2, (2) 

where M{H+5, |) is the magnetization at a field H+S (denoted by point A in Fig. 1(a)) from 
where the minor curve is initiated on the forward curve. Mml{H — 6, |) is the magnetization 
on the minor curve at a field H — 6, where it saturates as indicated by the point B in Fig. 
1(a). This procedure is based on the implicit assumption that the vortex state formed on 
the forward curve is an "equilibrium" state. This assumption is however inconsistent with 



the experimental observation by Wordenweber, Kes and Tsuei [T^, who showed that both 
current cycling and field cycling processes eventually establish a vortex state with a Jc higher 
than that on the forward curve. Such an observation indicates that the vortex state formed 
on the forward curve is metastable in nature. 



Tenya et al []T8[ have preferred a procedure given below, which is very similar to the one 



described above but using the minor curves of the type II described in Fig. 1(b): 



M,,{H) = [M{H -5,l) + Mml{H + S, T)]/2, (3) 

where H — 6 (point C in Fig. 1(b)) is the field from where the minor curve is initiated on the 
reverse curve and H + 6 (point D in Fig. 1(b)) is the field where it saturates. Mml{H + 6, | 
) is the saturated magnetization value on the minor curve. This procedure too has the 
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shortcoming similar to that in Eq. 2, viz., the vortex state on the reverse magnetization 
curve is actually a metastable state ||IO| , p!^ , p!6| . Moreover, not only are these recipes deficient, 
they also yield different conclusions, viz., an enhancement in equilibrium magnetization is 
observed in one case, whereas it is absent in the other. These ambiguities point to the 
need to evolve a more satisfactory procedure to arrive at a unique and stable vortex state 
unambiguously and determine the equilibrium magnetization assuming the stable state to 
be the equilibrium state. 



III. MODEL FOR HISTORY EFFECTS AND METASTABILITY 

Ravikumar et al incorporated the history dependence in the macroscopic critical 
current density Jc by postulating, 

J,{B + AB) = J,{B) + (|A5|/5,)(Jf - Jc). (4) 

where the critical current density Jc{B) is a macroscopic representation of a particular 
metastable configuration of the vortex lattice at a field B. Eq. 4 describes how the vortex 
state evolves from one metastable configuration to another. An important assumption of 
this model is the existence of a stable vortex state with a critical current density J**, which 
is unique for a given field and temperature. B^ is a macroscopic measure of metastability 
and describes how strongly could be history dependent. In the limit of Br tending to 
zero, however, this model reduces to the standard critical state model for which Jc (= J^*) 
is independent of the magnetic history. It can be seen from Eq. 4 that a metastable vortex 
state with 7^ Jc^i can be driven into a stable state by merely oscillating the field by a small 
amplitude (see Fig. 1 of the Ref. 20). In the PE region, the energy barriers between different 
metastable vortex configurations are much greater than the available thermal energy. The 
field cycling allows the vortices to move and explore the energy landscape and thereby 
rearrange in to a vortex configuration closer to the stable state. In the next section, we 
will demonstrate this experimentally and show that the stable state obtained is indeed 
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independent of the magnetic history. 

In the hmit AB —>■ 0, Eq. 4 can be rewritten in the form, 



±dJjdB = (Jf - J,) /Br, (5) 

where upper and lower signs are apphcable in the cases of increasing and decreasing local 
field B, respectively. In each case, the Jc(-B) can be obtained by solving Eq. 5, provided 
the functional form of Jl*{B) and Br{B) are known. We assume for J^^{B) and Br{B) the 
following forms used in Ref. pOf for calculating the minor magnetization curves: 



and 

Br{B) ^{B- i^oHio^rii^oHp - B)" for Hi^^ < 5//xo < 

~ otherwise (7), 

The first term in Eq. 6 is the field dependence of well below the peak and the second term 
reflects the peak in J^* vs B. Br{B) in Eq. 7 accounts for the observed history dependence 
in Jc in the PE region. 5^ = in the field ranges H < Hiow and H > Hp signifies that Jc 
is independent of the magnetic history and is always equal to the J^*. For the two limiting 
cases, H < Hiow and H > Hp, the intervortex interaction and the flux pinning are dominant 
respectively and therefore the stable state is readily accessed by the vortex lattice. The 
values of the different parameters used in this paper are listed in the caption of Fig. 2. 
Jc{H t) [Jc{H I)] is calculated by numerically solving Eq. 5 with the upper (lower) sign 
with the initial condition Jc{H t) [Jc{H |)] = Jf{H) at some field below Hiow (above Hp). 
In Fig. 2(a), we present an evaluation of Jc{H ]) and Jc{H [) which obey the inequality 
Jc{H t) < Jf {B) < Jc{H [). It was earlier interpreted that the vortex state formed on the 
decreasing field cycle is a supercooled disordered state . In other words, the vortex state 



formed in decreasing field (from above Hp) retains the memory of the vortex correlations 
from the previous fields. In analogy, we can argue that the vortex state formed on the 



increasing field cycle is a superheated ordered state. Both of these states are metastable in 
nature. As argued above, they can be driven into a stable state by oscillating the external 
field by a small amplitude. 

The magnetization hysteresis loop corresponding to Jc{H f) and Jc{H |) are shown in 
Fig. 2(b). Note the asymmetry in the hysteresis, usually observed in experiments. For 
a comparison, we also plot the magnetization hysteresis loop one would obtain within the 
framework of Bean's critical state model with Jc = J^* (applicable in the limit Br 0) 
which is symmetric in the forward and reverse field cycles, as shown by the dotted line in 
Fig. 2(b). Details of the magnetization calculation are described in Ref. 20. The minor 
magnetization curves of the types I and II calculated in the slab geometry, are shown in 
Fig. 3(a) and Fig. 3(b) respectively. They clearly mimic the behavior seen in experiments. 
We assumed Meq{H) = 0, in calculating these magnetization curves. We note that the 
calculated curves in Fig. 2 and Fig. 3 are not quantitative fits to experimental data, they 
only serve to illustrate the qualitative features of the observed data. 

In Fig. 3(c), we show M*^{H) determined from the calculated minor curves of the type 
I and type II following Eq. 2 and Eq. 3, respectively. The test of the self-consistency of 
these procedures lies in reproducing the original form {M^q = 0) assumed in the calculation. 
M*g{H) obtained from these two procedures are not only inconsistent with each other, but. 



also, do not conform to the original assumption that M^g = [|21|. The procedure of Eq. 
2 indeed produces a peak like structure in M*g{H) which has been shown earlier from an 
analysis of experimental data in 2H — NbSe2 following the same recipe . On the other 
hand, the use of Eq. 3 proposed by Tenya et al yields no variation in M*^ vs H across 
the PE region. Fig. 3(c) illustrates the unreliable and ambiguous nature of these recipes 
noted above and thus points to the need for a consistent approach in order to overcome their 
difficulties. 
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IV. EXPERIMENTAL RESULTS AND DISCUSSION 



In this section, we will show experimentally that repeated field cycling drives any 
metastable state into a stable state, which is unique at a given field We study the 

minor hysteresis loops traced by repeated field cycling and infer from these measurements 
the critical current density J^* and the equilibrium magnetization M^q of the stable state. 

DC magnetization measurements have been carried out using a Quantum Design (QD) 
Inc. SQUID magnetometer (Model MPMS5) in the peak effect region of a 2H-NbSe2 single 
crystal (T^ ~ 7.25 K) with the field applied parallel to its c-axis. The crystal is of approximate 
dimensions {a x b x c) 4mm x 5mm x 0.43mm. As stated earlier, the peak effect in Jc is 
manifested as the anomalous enhancement in the magnetization hysteresis (c.f. Fig. 1). 
The magnetization hysteresis has been studied at different temperatures from 6.7 to 6.95K. 
Magnetization hysteresis data at 6.95K was measured using a 2 cm full scan length, and the 
data at the other temperatures was obtained using the half-scan technique [|12|,^ to avoid 
artefacts arising due to field inhomogeneity experienced by the sample along the scan length. 
In the temperature range investigated, Jc at the peak field Hp decreases with decreasing 
temperature (see Table 1). 

Fig. 4(a) depicts a part of the hysteresis loop at 6.95K, constituting M vs H curves 
in the increasing (forward) and decreasing (reverse) field cycles measured with a 30 sec 
wait time at each field. We identify the onset field H^i of the PE on the forward curve, 
where M begins to decrease sharply. The field Hp marks the field at which magnetization 
hysteresis is maximum. In Fig. 4(a), we show the points A, B, C and D from where the 
minor hysteresis loops are initiated. A(C) and B(D) are at a field H < H^i {H > H^i) 
on the forward and reverse curves, respectively. Minor hysteresis loops starting from both 
forward and reverse curves are recorded at different fields (spanning the peak region) by 
repeatedly cycling the field by a small amplitude AH. The interval AH is chosen such 
that it is above the threshold field required to reverse the direction of the shielding currents 
throughout the sample. From the critical state model, we understand that magnetization 
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values must always remain confined within the forward and reverse magnetization curves, 
which constitute the so called magnetization envelope. Further, the M — H loop in each 
field cycle must retrace itself. 

In Fig. 4(b), we show the minor hysteresis loops (MHLs) measured by repeatedly cycling 
the field, starting at point A (H < H^i) on the forward curve. These MHLs in different 
field cycles retrace each other indicating that the Jc does not change with field cycling. 
Therefore, we conclude that the vortex state is in a stable configuration. In contrast, the 
MHLs shown in Fig. 4(c) (continuous line with data points omitted) starting at B (if < i?^) 
on the reverse curve, show shrinkage effects, with each successive field cycle and finally MHL 
collapses into the minor loop started from point A (open circles) which is replotted in Fig. 
4(c). This suggests that the vortex configuration at point B is metastable with a Jc > Jf. 
Repeated field cycling causes the Jc to fall towards the stable stable value as reflected in 
the reduction of the width of the MHL with each successive field cycle. It is remarkable 
that the minor loops starting from both A and B merge into precisely the same loop within 
the experimental accuracy. This clearly reaffirms the basic assumption of the model that 
there exists a unique stable state with critical current density Jf, independent of the initial 
vortex state from which it evolves. 

We now focus on the behavior of MHLs which start from a field H > H^i. As shown 
in Fig. 4(d), the behavior of the minor loops starting at point C is quite different from 
those started at point A. The increasing field leg of the MHL moves away from the forward 
magnetization curve in the first field cycle itself and remains outside the magnetization 
envelope for subsequent field cycles. This clearly suggests that, for H > H^i, the vortex 
configuration even on the forward magnetization curve is metastable. However, the behavior 
of the MHLs starting at point D on the reverse magnetization curve is very similar to the 
behavior of those that start at point B, i,e., the MHL shrinks with each successive field 
cycle (continuous line in Fig. 4(e)). The data in Fig. 4(d) is replotted in Fig. 4(e) (open 
circles connected by dotted line), which suggests that the MHLs starting from both C and 
D collapse into the same final loop (MHL). Firstly, the data in Fig. 4 clearly suggest the 
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metastable nature of the vortex configuration for fields above H^i both on the forward and 
the reverse magnetization curves. Further, the eventual MHL obtained on repeated field 
cycling is independent of the initial vortex configuration. We note that the metastable state 
on the forward magnetization curve settles into the stable state much faster than that on 
the reverse curve. This might imply that the vortex configuration on the increasing field 
cycle is closer to the equilibrium configuration. 

The data in Fig. 4 yield the following inequalities for the critical currents in the different 
field ranges: (i) For H < H^i, the vortex configuration is stable in the increasing field cycle 
while at the same field value, it is highly metastable in the decreasing field cycle. This can 
be summarized by the inequality JdH t) = Jc^{H) < Jc{H |); (ii) For H^^ < H < Hp, the 
vortex configurations in both increasing and decreasing field cycles are metastable, with the 
critical currents obeying the inequality, Jc{H t) < Jc^{H) < Jc{H |); (iii) For H > Hp, 
Jc{H ]) = Jc{H I) = Jf{H). These observations are in accordance with the model pO[ (cf. 



Fig. 2(a)). We thus assert that the Eq. 2, proposed by Roy and Chaddah |T^ is applicable 
only for H < H^^. It is unsatisfactory for H^^ < H < Hp, as the vortex lattice on the forward 
curve is in a superheated vortex configuration which is more ordered (but metastable) than 



the stable configuration. Eq. 3, as proposed by Tenya et al |TB| is not appropriate in any 
of the field ranges because the vortex states produced on the reverse curve are supercooled 
vortex configurations [p!4| , ^ which are more disordered than the corresponding stable states. 

Fig. 5 shows the M-H loop at 6.9K constituting the forward and reverse magnetization 
curves (dark line with data points omitted) indicating H^i and Hp. Note the asymmetry 
(also seen at 6.95K) in the forward and reverse magnetization curves which is the hall mark 
of the peak effect. We also measured the MHLs by repeatedly cycling the field starting at 
different points on the forward and reverse curves. The saturated MHLs are again found to 
be independent of the initial vortex state just as for 6.95K. The locus of magnetization values 
on the increasing and decreasing field legs of the saturated MHLs measured at different fields 
are also plotted in Fig. 5 (open circles connected by dotted line). This observed behavior 
is in excellent qualitative agreement with that expected from the model in Ref. 20 (see 
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Fig. 2(b)). The locus of saturated magnetization values corresponds to the "table" or the 
"quilibrium" vortex configuration at different fields. 

Having established the existence of a history independent stable state we determine the 
critical current density J^* and the equilibrium magnetization M^q state at each field from 
the saturated MHL J^* and Meg are given by, 



Jl\H) = {Mst{H i) - Mst{H ])]/2g^^^R, (8a) 

M,q{H) = [M,t{H t) + M,t{H i)]/2, (86) 

where Mst{H ]) and Mst{H I) are the magnetization values on the increasing and decreasing 
field legs of the saturated MHL. J^* vs H and M^q vs H data at 6.95K are plotted in Fig. 6(a) 
and Fig. 6(b), respectively. M^g exhibits a sharp increase between H^i and Hp signifying 
an increase in the equilibrium fiux density. This is reminiscent of the characteristic of M^q 
across the FLL melting transition observed in cuprate superconductors p5| , |26| . We argue 



that the change in M^q indicates a first order transition in the FLL from an ordered solid to a 
pinned amorphous state [|19| presumably analogous to a Bragg Glass to Vortex Glass/pinned 
liquid phase transition The increase in M^q coincides with the increase in J^* near the 



onset of the peak effect and spans the field range between H^i and Hp. In Fig. 7(a) and 



7(b), we present the M^q vs H and Jf vs H data respectively, at 6.9K. We note that the 

igq i^OrrelateS with a onaip ihl-iccidc m uctwccn J^-ipi aiiu. -i-ip. 



sharp change in M^q correlates with a sharp increase in J^* between if and Hp. We also 



present the AMgq values obtained at different temperatures in Table 1. 

It is important to understand the nature of the vortex state in the transition region 
Hpi < H < Hp. One of the well known pictures is the collective pinning scenario 0, where 
the loss of long range order is expected to permeate uniformly throughout the sample. On 
the other hand Paltiel et al |^ have recently proposed a picture where the disordered phase 
enters through surface imperfections and coexists near the surface with the ordered phase 
of the bulk. They argue that the boundary between the disordered region and the ordered 
region moves into the sample as the temperature (or field) is increased towards Tp (or Hp). 
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Further possibility is the coexistence of ordered and disordered phases, with an intricate 
geometrical connectivity of these phases. Irrespective of the particular picture used, our 
experiments demonstrate a specific and an unambiguous procedure, viz., subjecting the 
sample to a field cycling, to produce a unique stable state (in a macroscopic sense) across 
the peak effect region. 

We consider this stable state as a pinned equilibrium state, and estimate equilibrium 
magnetization and the free energy difference or entropy change when the vortex lattice 
changes from an ordered to an amorphous state. As per the Clausius-Clapeyron relation 
26| , |29[| , the entropy change per vortex per inter-layer distance 4 Angstroms) in the 



2H - NbSe2 system [19 



As = -{AM,jH,){dH+/dT){(j)od/kB), 

where dH^JdT ~ dHp/dT ^ —0.65 T/K. The value of As estimated at different tempera- 
tures is tabulated in Table 1. Incidentally these values are comparable to the entropy change 
reported across the FLL melting transition in high Tc cuprates. 

An important question that can arise is whether the entropy change can be observed in 
thermal measurements such as specific heat vs temperature. We recall that the metastability 
in the vortex state is much greater in temperature scans in a fixed magnetic field |]14| . 



Repeated cycling of the field by a small amplitude may be necessary to produce the "stable" 
or "equilibrium" state before a thermal measurement is carried out at each temperature. 



V. CONCLUSIONS 

In this paper, we have presented a study of the different metastable vortex configurations 
occuring in the peak effect region of a weakly pinned superconductor 2H — NbSe2 through 
magnetization measurements. Each metastable vortex configuration is characterized by a 
critical current density Jc which is strongly dependent on the magnetic history. It is also 
shown that any metastable vortex configuration obtained under given field historys can be 
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driven into a stable configuration by repeated field cycling. This stable configuration has 
a critical current density J^*, uniquely determined by field and temperature as postulated 
in a recent model [Q. Field cycling appears to act as an effective temperature to drive 
a metastable state into the stable state, even when thermal energy itself is inadequate to 
sample the phase space and access the stable state. 

The method of recording minor hysteresis loops described here allows us to determine the 
pinning and equilibrium properties of the stable vortex state satisfactorily. Our equilibrium 
magnetization data clearly suggest that the transition of the vortex lattice from an ordered 
state to a disordered state is first order in nature. The smearing of the transition, i.e., 
the width of the transition region may be a manifestation of the spatially inhomogeneous 
pinning of the system. The J^* data suggests that the loss of quasi-long range order in the 
vortex lattice also spans the same field window as the magnetization jump. In the collective 
pinning picture, this amounts to correlation volume of the vortex phase decreasing in this 
regime and the FLL becoming completely disordered above Hp or Tp. The precise coincidence 
of the Jc anomaly with the equilibrium magnetization anomaly further illustrates the self 
consistency of the procedure developed here. It would be interesting to compare the nature 
of this disorder-driven transition in systems with different types of pinning, e.g. high density 
of point pins versus low density of extended pins to further understand the nature of this 
presumably disorder induced phase transformation. 

The authors thank Dr. K. V. Bhagwat, Dr. T. V. Chandrasekhar Rao, Dr. P. K. Mishra 
and Mr. M. R. Singh for discussions. 
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TABLES 



TABLE I. Superconducting parameters in 2H — NbSe2 

T(K) Hp{mT) MoAMeg(;ur)/47r MHp){A/m^) Asiks) 

6.95 105 3.8±0.4 54 x 10^ 13.6±1.4 

6.90 136 5.0±0.4 36 x 10^ 13.8±1.1 

6.85 170 1.3±0.4 26 x 10^ 2.9±0.9 

6.80 202 1.9±0.4 17 x 10^ 3.5±0.7 
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FIGURE CAPTIONS 



Fig.l. : Typical magnetization hysteresis loop observed in the peak effect region of a su- 
perconducting 2H — NbSe2- In the panel (a), the minor curve obtained by decreasing the 
field from the point A, corresponding to a field {H + 5) on the forward magnetization curve 
is shown to saturate at the point B, which corresponds to the field {H — S). Magnetization 
values at A and B are M{H + 5, t) and Mml{H — 5, 1) respectively (see text). In the panel 
(b) the minor curve obtained by increasing the field from the point C, corresponding to a 
field {H — S) on the reverse magnetization curve, saturates at D (if + 5) and corresponds 
to a magnetization value Mml{H + 5, t). 

Fig. 2: (a) Calculated critical current densities Jc{H f) and Jc{H I) in the increasing 
and decreasing field cases, respectively. These are compared with the stable critical current 
density (dotted line). In this calculation, we have used Hiow — 0.05T, Hp = O.IT, Jd 
= lO^A/m^, = 20Jci, Hi = 0.12r and Hw = 0.008T [20]. (b) Magnetization hysteresis 
loop corresponding to the Jc values shown in (a). The hysteresis loop that would be obtained 
within the framework of critical state model, i.e., in the limit oi B^. ^ is also shown in the 
in this panel as dotted line. The inset shows the functional form of Br which is non-zero in 
the field range Hiow < H < Hp. 

Fig. 3: Calculated minor curves of type 1 and type 11 are shown in pnaels (a) and (b), 
respectively. In the panel (c), we show the M*^ vs H obtained using Eq. 2 and Eq. 3, 
respectively along with the original form, M^g — 0, assumed in the calculation of the minor 
curves. 

Fig. 4: (a) A part of the magnetization loop (forward and reverse curves) measured at 
6.95K on a 2H — NbSe2 single crystal. Also indicated are the characteristic fields, H^i and 
Hp. We indicate A and B {H < H^i) and C and D {H^i < H < Hp) starting from which 
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the minor hysteresis loops are measured, (b) Minor hysteresis loops started from point A 
(open circles). In different field cycles, they are seen to retrace the same loop, (c) The MHL 
started from B (continuous line) shrinks with each successive field cycle. The increasing and 
decreasing field legs of the first and second cycles are numbered. After five field cycles, the 
hysteresis loop is seen to merge with the loop shown in (b), which is replotted (open circles), 
(d) Minor hysteresis loops started from point C (open circles). In the first field cycle itself, 
increasing field leg of the MHL moves away from the forward curve and remains outside the 
magnetization envelope for the subsequent field cycles, (e) The minor loops starting from 
D (continuous line) are seen to collapse onto the loop shown in (d), which is replotted. 

Fig. 5: Magnetization hysteresis loop of 2H — NbSe2, recorded using half scan technique [22] 
at 6.9K (continuous line). The open circles are the saturated magnetization values obtained 
after repeated field cycling. H^i and Hp are also marked. The locus of the saturated 
magnetization values is shown as a dotted line. 

Fig. 6: (a) Stable critical current density J^* in the field range 80 mT < /iqH < 105 mT. In 
the inset, we show the Jf vs H in the entire field range. Filled triangles and open circles 
correspond to the values obtained from the MHLs intiated from the forward and reverse 
magnetization curves respectively, (b) Equilibrium magnetization M^g as a function of field 
at 6.95K. Note that the sharp change in Mgg coincides with the PE onset field if^. 

Fig. 7: (a) Critical current density vs H and (b) Meg vs H data obtained at 6.9K. 
Note that the smeared jump in Mgg vs H, as marked by the double sided arrow, agrees 
precisely with a smeared jump in critical current density in the peak regime. See text 
for a discussion. 
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